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Primary 




science for GBM is detection of Gamma-ray Bursts 


CTIME 

CSPEC 


- 0.512 s time resolution, 8 channels 

- 4.096 s time resolution, 128 channels 


Launched June 1 1 , 2008 














Large Area Telescope 
(LAT) 

00 GeV 

Gamma-ray Burst Monitor 

12 Nal detector. 

8 keV - 1000 keV 
126 cm 2 , 1.27 cm 

Triggering, localization, spectroscopy. 

2 BGO detector. 

200 keV - 40 MeV 
126 cm 2 , 12.7 cm 
Spectroscopy 

Bridges gap between Nal and LAT. 






GBM Earth Occultation Project 

PI Colleen Wilson-Hodge 

All sky X-ray monitor of known sources from 8 

keV - 1000 keV 


Source Database 
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Conceptually simple 




Fermi 


Gamma-ray Burst Monitor 


GBM Earth Occultation Method 


Predict occultation times 


In practice.... 


Determine detectors 
viewing source of interest 

Fit to each detector and 
energy channel 

- Background model 

- Model count rates for 
each source 

• Detector responses 

• Assumed energy 
spectrum 

• Atmospheric 
transmission 

Compute best scale factor 
for all detectors to 
estimate fluxes. 


650 r 


600 


550 


500 


450 - 


400 - 


350 1 


Fit Window 

(240 s) 
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Occultation Time 


The time where the probability that a 100 keV gamma ray from the 
source will pass through the atmospheric column is 50% 


Atmospheric Transmission function 

T(E ph ,t) = exp[- fi(E ph A(h(t))] 


A(h(t )) 


mass attenuation coefficient of gamma rays at photon energy^,}, in air 

air mass along the line of sight at a given altitude hit ) based on the U.S. Standard 
Atmosphere (1976) 



• Each detector which views the source of interest within 60 degrees of the detector 
normal is included in the fit 

• Observed count rate model for each detector is: 

n 

r(t, E ch ) = b 0 (E ch ) + h(E ch ) * (t - t 0 ) + b- 2 {E ch ) *{t- 1 0 )' 2 + ^ a.i(E ch ) * Si(t, E ch ) 

i=l 

bo(E c h), bi(E ch ), 62 {E c h) = Quadratic background coefficients 
a % {E c h) — Fitted scale factors for each source model 

Si(t,E ch ) —Source models for source of interest and all other sources included in the fit window 

S(t , E c h) — R(E p hi E c h, t) f)) * f(E p h)dE p t^j 

f(E ph ) — Assumed source spectrum 

T(B ph , t) = Atmospheric transmission 

R(E p h, E c h, f) =Time dependent detector response 




Fitting 





Flux Measurements 

Each energy channel and each detector is fitted independently 

F(E c h) =a(E c h) * f f(E p h)dE p h 

J Eph 

a(Ech) = Weighted mean of scale factors for each detector 



Interfering sources in fit window 

Each source in the database is identified as: 


> 500 mCrabs 
>150 mCrabs 

> 50 mCrabs 


• Strong - Always include in fit out to 90 degrees 

• Moderate - Always include in fit out to 60 degrees 

• Weak - Always include in fit out to 40 degrees 

• Quiescent - Never include in fit unless it is flaring 


Flare database 

Public SwiftfB AT transient monitor data 

• 50 mCrabs < Source <150 mCrabs - Weak 

• 150 mCrabs < Source < 500 mCrabs - Moderate 

• Source > 500 mCrabs - Strong 

If an interfering source meets the criteria for any detector it is included for all detectors 


Additional Considerations 



Eclipsing sources 

10 sources in the catalog are eclipsing 

Sun-Solar flare database 

• Class M or X flares - Strong 

• Class C flares - Moderate 

• Class B flares - Weak 

Pre-Filtering data 


Usually removes class M and X flares as well as SAA entrances and exits 



Systematic Effects 

• Accuracy of assumed source spectral model 

• Heavily tested and researched; -3 power-law 

• Large variation in background 

• Pre-filtering of data 

• Duration of the occultation transition 

• High latitude sources; Limited to 20 seconds 

• Inaccuracies in the detector response matrices 

• Remove steps for all possible solar panel blockages 

• Occultation limb geometry 

•52 day precession; Flare database 

• Nearby sources 

• Exclude steps if bright source is within 8 s of occultation 
time 
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Gamma-ray Burst Monitor 


Post-Filtering 


Occultation steps are removed if: 

The source of interest occults within 8 s of a bright source 
The occultation lasts for longer than 20 s (high latitude sources) 
The space craft is rapidly slewing with a spin rate> 0.004 rad s 

Individual steps are > 10cr or > 3.5<r from the mean if sources 
intensities reach 150-500 mCrab or < 150 mCrabs respectively 
The time of the fit window is associated with a solar flare 
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Ghost Source Analysis Systematic Errors 
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^ ^ &total & stat 

o total = Width of flux distribution 
k = l-O/asig Scale factor 

® sys ® total ® stat 


Systematic Error Estimates for GBM Earth Occultaiion Analysis 
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Sco X-l 


Comparison Between GBM and Swift ! BAT 


GBM 12-50 keV 


4U 1700-377 


A 0535-262' 




CRAB 


Cyg X-3 


GSR 1915+105 •> j 

44iiMi8i4.il 




Vela X-l 


Swift/BAT 15-50 keV 


2-4 day averages 



Three Year Fermi/G BM Earth Occultation Catalog 

• Source Name 

• Ra & Dec 

• Category (A, B , T, P, N, I) 

• 3 Year Average Flux (mCrabs) 

- 12-25 keV 

- 25-50 keV 

- 50-100 keV 

- 100-300 keV 

• Significance 

- 12-50 keV 

- 12-300 keV 
•Type 

The Astrophysical Journal Supplement Series, 201:33 (22pp), 2012 
August 






Detection Criteria 

• Significance exceeds 5 or 3 .5 sigma (Category A and B respectively) 

• Detected in the transient search (T) at 5 or 3 .5 sigma 

• Detected in the orbit folding search at 5 or 3.5 sigma (P) 

Non-Detections 

• Significance less than 3.5 sigma (Category N) 

• Significance is negative (Category I) - only 6 sources 



Summary of Results 


3 Year Catalog 

209 Sources (99 detected- A) 

40 LMXB/NS 
31 HMXB/NS 
12BHC 
12AGN 
1 Star (Sun) 

1 TDE (SWIFT J 1 64449 .3+57345) 
1 Pulsar/P WN (Crab) 

1 Galaxy Cluster (Oph Cluster) 


Current 

215 Sources (104 detected- A) 

40 LMXB/NS 
34 HMXB/NS 
14BHC 
12AGN 
1 Star (Sun) 

1 TDE (SWIFT J 1 64449 .3+57345) 
1 Pulsar/P WN (Crab) 

1 Galaxy Cluster (Coma Cluster) 





Transients Seen with Earth Occultation 
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Sources Detected Above 100 and 300 keV 






I 


0-0* 

0-05 

0.04 

0.03 

0.02 

15.01 


i 00-300 k*V 





300-500 k*v 


53000 55200 

Tim? (MJD} 


55400 


5*000 






5*600 


Crab 


1 00 -300 ktv 


300-300 keV 


550150 55200 

Time {MJD] 


Centaurus A 


Cygnus X-l 


5 00-300 kev 


55400 




O D . 




Cd" 0 ' 2 


U-Q.* 


X 2 r 


300-500 keY 




b*BVQ 


£5000 ££200 

Throe (MJD) 


S6*M 


Case et al. 


IE 1740-29 






Flux (Crab) Flux (Crab) 



Sources Detected Above 100 and 300 keV 
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Crab Flux Decline 
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Crab Spectrum 


1 — j - 
r t =2.04 
62 keV 
= 2.18 


(Jourdain and Roques 2009) 
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More spectral analysis coming to our web site in the near future 


Flux (Crob units) 



G. L. Case et al. 
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Monitoring of Cygnus X-1 During the 2010-12 State Transitions 

with the Fermi GBM 
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Brief History of X-ray Occultation 



^^2008- Earth Occultation measurements with GBM 




Occultation Coverage 


The diameter of the Earth seen from Fermi is ~ 140°, so roughly 30% of 
the sky is occulted by the Earth at any one time. The precession of the 
orbit, caused by the oblateness of the Earth means that the entire sky 
is occulted every ~53 days. 


Fermi Orbital 
Information 

•Scans the sky pointing 
+35° one orbit then slews 
pointing -35° 

•Rolls about z-axis to keep 
solar panels facing sun 



V, 



| 3 , the ‘elevation angle’, is measured at the geo-centre between a 
source crossing the Earth’s limb and the orbital plane. 




Imaging with Earth Occultation 


• Tomographic method that builds images from projections of 
Earth’ s limb (edge) on the sky 

• Detectors uncollimated so detector count rate measure of sky 
brightness 

• For detector facing Earth, time derivative of count rate is 
proportional to line integral of flux along Earth’ s limb 

• During a precession period, elevation angle, 0(t), changes so 
flux is measured at different angles w.r.t Earth’ s limb giving 
different projection angles 

• All available projection angles sampled over course of 
precession period 
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Limbs for CRAB on MJD 54690 = 16.6 



Limbs for CRAB on MJD 54700 = 45.1 



Limbs for CRAB on MJD 54720 = 13,7 



Limbs for CRAB on MJD 54730 0 * = -3,3 



Limbs for CRAB on MJD 54710 = 41.0 





Limbs for CRAB on MJD 54690.0 to 54690. 
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Occ. Limbs for 3C273 on MJD 54690.0 to 54690.0 
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Occ. Limbs for NGC4151 on MJD 54690.0 to 54690. 
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Motivation for Imaging 


• Look at high energies to settle discrepancy between 
EOT vs. EBOP 

• Construct complete catalog to reduce a cause of 
systematic errors 

- Flux from unaccounted-for sources attributed catalog sources 

- Example with 1 E 1740-29 and transient GX 1+4 during flare 
MJD 55708-55774 
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Imaging with IDEOM 


• Grid of virtual sources at 0.25° 
spacing, -660,000 sources 

• Occ. times predicted 

• 4-min window of count-rate data 
centered on occ. time 

• Select detectors view steps < 75° 

• Detector response correction to 
account for angle between source 
and detector normal and energy 
efficiency 

• Convert rise occ. steps to set 
steps; rotate about center 

• Average together data from all 
detectors for all windows for a 
given potential source 



Image: Meegan et al. 2009 


Imaging Technique 


• Apply differential filter to 
average window 

• f a = inner boundary = 3 

• f b = outer boundary = 8 

• Sums 8 bins of data, skips 6 
bins, sums the next 8 bins, 
and then takes the 
difference between the 8-bin 
totals 

• Fit polynomial to center 
region 

• Fit a spline function to the 
side areas for background 

• Take the difference 
between the polynomial and 
spline at the occultation time 
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Implementation 


-660,000 sources 

Parallel machines to process multiple sources 
at once 

Use LSU HPC Tezpur and LONI Queen Bee 
machines 

Split catalog into 3 roughly equal sections 

Takes -6 hours with 80 nodes (320/640 
processors) to run 1 catalog for a single 
energy band 




IDEOM Results 


28 precession periods (~ 4 years) for: 

- 1 2-50 keV 

- 50-100 keV 

- 1 00-300 keV 

Single precession period images for 3 energy ranges 
Images summed weighted by number of occ. steps 
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12-50 keV ~4 Years 
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Finding Sources with IDEOM 


Early searches cross-correlated features in images 
with BAT, INTEGRAL, and LAT 

- 16 sources add to catalog this way 

Systematic search 

- For virtual sources >3.5 a, 5°x 5° regions fit to 2D Gaussian 

• Centroid taken as source position 

• Peak value as source significance 

• Integrated significance from fit Gaussian widths 

- Candidate sources product of peak significance and 
integrated over stat. error > 10 

- Candidate sources analyzed with EOT 




Comparison to GBM Catalog 


12-50 keV: 

- 67 sources detected by GBM-EOT; 43 by IDEOM 

• 1 1 blended with brighter 

• 12 no feature 

• 2 with test statistic <10 

• 1 is A0535+262 

- Sensitivity limit 25 - 50 mCrab 

- Position accuracy ~ 0.4° 




Comparison to GBM Catalog 


50-100 keV: 

- 25 sources detected by GBM-EOT; 23 by IDEOM 

• 2 blended, 3 no feature, 3 IDEOM, but not sig. in GBM 

- GRS 1 758-258 and 4U 1 81 2-1 2 show flux out to 1 00 
keV so consistent with INTEGRAL 

- NGC 4388 flux consistent between IDEOM and 
GBM-EOT 

• Sensitivity limit ~ 22 mCrab 

• Position accuracy - 0.47° 




Comparison to GBM Catalog 

1 00-300 keV: 

- 6 sources by GBM-EOT; 7 by IDEOM 

• 6 GBM-EOT + GRS1 758-258 

• Sensitivity -32-54 mCrab 

• Position accuracy - 0.25° 

• Only GBM sources down to -20 mCrab 

- No unaccounted-for sources suggests 
EBOP background model inaccurate 




GBM-LAT Sources 


GBM and LAT give energy coverage 10 keV to 
300 GeV 

4 persistent sources detected by GBM and LAT 

- NGC 1275 

- 3C 273 

- Cen A 

- Crab 

Data from GBM, LAT, and COMPTEL 




Where do the “hard” photons 

come from? 


“soft” photons 
synchrotron emission 
from AGN jet 

“Hard” photons from 

- Inverse Compton 
scattering (1C) 

• Synchrotron self- 
Compton (SSC) 

• External inverse 
Compton (EC) 

- Proton Induced Cascade 
(PIC) 



1 



NGC 1275 


• Blazar at ~75 Mpc (z=0.01 7) 

• Detected in hard X-rays by 
HEAO 1 (1977-1979) 

• Not seen by CGRO (1991-2000) 

• Monthly variations seen by LAT 

• Best power-law spectral fit 
r=2.13 

• Best fit cut-off power-law l~=2.07 
and E c =42.2 GeV 

• LAT power-law plotted 

• Fits GBM data reasonably well 

• Flux due to SSC (Abdo et al. 
2009) 
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3C 


• Quasar at -700 Mpc (z=0.158) 

• Spectral hardening >1 MeV 

• Von Montigny et al. (1997) fit 
using SSC, EC, and PIC 
models 

- Best fit by PIC 

- Fit X-ray / y-ray data 
empirical model 



- Plotted with n,=0.7, r 2 =0.8 
and E b = 2.36 MeV 

- Consistent with data to -10 
GeV 


273 





Cen A 


Radio galaxy at ~4 Mpc 

Seen from radio to TeV energies; 
UHE cosmic rays 

CGRO showed spectral variability of 
increasing break E as source 
brightens 

LAT fit radio to GeV y-rays with SSC 

Reynoso et al. (201 1 ) fit radio 
through TeV with model using e- 
synch, p-synch, and 1C 

Reynoso et al. (201 1 ) overplotted 
Consistent with data out to ~10 GeV 



Crab 


Pulsar/PWN at ~2 kpc 

Emission to - 100 MeV from 
synchrotron 

E > 100 MeV from 1C, SSC 

Dashed line from INTEGRAL 
with r 1 =2.07, r 2 =2.23 and 
break at 62 keV 

Solid line from GBM with 
r 1 =2.06, r 2 =2.36 and break 
at 98 keV 

GBM fit in very good 
agreement with LAT data 
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Transient Search 


Wide F.O.V. for frequent monitoring; can observe source during 
entire flare 

1-day 3 o sensitivity in 12-25 keV -150 mCrab 
Use sliding average to increase sensitivity in search 
- 5-day, 9-day, and 19-day binning 

Searched 12-50 keV, 50-100 keV, 100-300 keV, 300-500 keV 
bands 

Joined consecutive bins >3o as 1 event 
Events with significance >7o 

Computed significance above median flux; keeping >6o 




Search Results 


Searched 4 years light curves for all GBM sources 
Detected 168 transient events from 65 sources 

- 7 above 50 keV 

- 1 above 100 keV; XTE J 1752-223 

Comparison to GBM flare database 

- Limited to database events at least 5 days 

- 280 events from 44 sources 

• Difference mainly due to periodic sources: SMC X-1 , LMC X-4, 
HerX-1, and EXO 2030+375; > 100 missed detections 

• 72 just from SMC X-1 and LMC X-4; at low declination 

- 10 sources in database not in search; only 1 >6o (X Per) 




GK 


• GBM flare avg flux ~51 mCrab; 
one of lowest detected from 
flare search 

• Intermediate polar system 

• Cooled mainly by thermal 
bremsstrahlung radiation 

• Spectral blackbody fit give kT e 

• Able to estimate lower limit 
mass 
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GX 339-4 


• Historically detected > 100 keV 

• Including marginal detection, 
GBM observations of 2 nd flare 
in good agreement with 
INTEGRAL observations 

• In hard spectral state 


F(E) = NE~ r e~ E,Bc 


• r=1.69 and E c =328 keV 

• Emission mostly from 1C by 
soft X-rays from corona 

• OR from 1C from base of jet 

• Jet and corona models diverge 
> 1 00 keV 




XTEJ 1752-223 


• BHB discovered Oct. 23, 2009 

• Due to solar constraints, GBM 
only instrument to cover 
duration of hard state 

- RXTE missed 2 months 

- INTEGRAL ~ 1 month 

• RXTE observations shortly 
after outburst show hard state 


12-50 keV 
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50-100 keV 
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- r=1.24 and E c =133 keV 

• No available jet/corona models 
to compare 


100-500 keV 
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Conclusion 

EOT applied to GBM to monitor in hard X-ray/soft y-rays 
All-sky imaging with IDEOM for 

- ~4 years of data in 3 energy bands 

- ~45 sources detected 

- 16 added to GBM catalog added by cross-correlating with other catalogs 

- Position accuracy -0.5° 

- No additional sources detected > 100 keV so likely EBOP background 
inaccurate 

Spectral analysis of GBM-LAT sources (NGC 1275, 3C 273, Cen A, Crab) 
show data consistent with models from lit. 

Transient search 

- In 4 years 168 events from 65 sources; 7 >50 keV; 1 > 100 keV (XTE 
J 1752-223) 

- Able to detect flares down to ~ 50 mCrab 

- No solar constraints enable important observations (XTE J 1752-223) 

- High energy capability help to distinguish between jet and corona BHB models 




Systematic Errors 


• Occultation Limb Geometry 

- Occ. time corresponds to arc 
along Earth’ s limb 

- If unmodeled source occults at 
the same time can cause 
inaccurate fluxes 

- Projection of Earth’ s limb on 
sky repeats ~ 53 days creating 
periodic errors 

- Add newly discovered 
transients to catalog (BAT/ 
MAXI) 

- Imaging with IDEOM to look for 
sources 

• Nearby Sources 

- Unable to properly assign flux 
when known sources occult 
within ~8 sec of each other 
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Previous Imaging with Earth 

Occultation 


Occultation transform imaging 
(Zhang et al. 1993,94,95) 

- Reconstructs projection data 
using inverse linear Radon 
transform; used in CT scans 
and medical tomography 

- Uses Maximum Entropy 
Method (MEM) to locate 
sources 

- Able to image 20-300 keV 

- 30° x 30° regions along GP 
looking for transients 

- Limited to small regions of sky 
where limb approx, linear 

- Doesn’t account for bright 
sources outside f.o.v 


• Likelihood Imaging Method for 
BATSE Occultation data 
(LIMBO) (Shaw et al. 2004) 

- Predetermined grid of positions 

- Background subtract data 

- Pass through differential filter 

- Calculate likelihood flux due to 
a location compared to null 
result 

- Able to image entire sky, ~ 500 
days in 25-160 keV 

- Limited by 2° separation of grid 

- Only single energy range 





http;//heastro .phys .lsu .edu/gbm/ 




